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Preface 

Like in the laser beam, where the Bose-Einstein statistics of photons results in its high 
momentum coherence, in the production of large number of pions this statistics may lead to 
strongly enhanced emission of pions with the same charge. Because of that, in production of a 
hundred pions, the events without neutral pions or, on the contrary, without the charged ones, 
extremely rare from point of view of classical statistics, obtain appreciable probability and can 
be observed. They were observed in the cosmic rays experiments (Centauro events) at energy 
~ 1000 TeV. This energy region will be covered and studied in detail in the LHC experiments. 

The consequences of the Bose-Einstein statistics for pion creation and their charge dis- 
tributions is being discussed in the literature for a rather long period of time (see, e.g., the 
paper [I\ in the list below). A possibility of formation of a large domain of disoriented chiral 
condensate was indicated in the work [//] and then in [III,IV\. The more detailed lists of refer- 
ences devoted to this subject can be found in the more recent publications, see, e.g., the papers 
[V,VI,VII\. 

As far as we aware, our paper [a] was one of the first where the pion charge distributions 
were calculated, proceeding directly from the Bose-Einstein statistics, and where the analogy 
with the induced photon emission was strongly emphasized. However because of some technical 
reasons this paper was never sent to any physical journal. Its short version was published in 
[6] (in Russian) and is practically inavailable. Therefore we put now the paper [a] in the LANL 
archive. No changes have been done in this text except for correcting some misprints. 



[a] V.A. Karmanov and A.E. Kudryavtsev, The Centauro events as a result of induced pions 
emission, Preprint ITEP-88, Moscow, 1983. 

[h] V.A. Karmanov and A.E. Kudryavtsev, The charge distributions of creation of coherent 
pions, In: "Nucleon-nucleon and hadron-nucleus interactions at intermediate energies", Proc. 
of third symp. of LINP, Leningrad, 1986, p. 558 (in Russian). 
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Abstract 



It is shown that only in the case of n pion production amphtude (at n > 10), which is 
symmetric enough relative to the momenta permutations, there is appreciable probability 
of events with nothing but charged or nothing but neutral pions. Significant enhancement 
of the probability of these events in comparison with that follows from Poisson distribution 
is caused by induced pions emission, in complete analogy with the phenomenon of the 
induced photon emission in QED. 

1 Introduction 

Overwhelming majority of theoretical pictures describing the high energy processes at high 
multiplicity leads to Poisson distribution for probability of emission a given number of pions, 
provided the average number of pions is fixed. In more sophisticated models this probability 
is given by the sum of Poisson distributions - see papers and references therein. In all 
the distributions of that kind the average number of neutral no and charged pions rich are 
determined by the following relatons: 



Notice that the experimanteal data obtained by the existing accelerators {^/s < 540 GeV) do 
agree with that theoretical distributions. On the other hand, in cosmic rays the events were 
found with large multiplicity (n ~ 10-i- 100), in which the neutral pions are absent p|, ^. Their 
energy is about Eiab ~ (1 -h 2)10^ TeV. These events were called Centauros. Below we shall 
suppose that charged particles in Centauro events are the pions. 

Notice that it is practically impossible to interpret Centauro events as a fluctuations of 
Poisson distribution - the probability w{7c~^7t~) for such a fluctuation is too small. For example, 
for n = 100 we get: 



As it will be shown below, the charge distribution in a system of the given total num- 
ber of pions n is closesly connected with symmetry properties of the production amplitude 
M(pi, Ti, . . . ,pn, Tn) relative to permutations of isotopic indices ri, . . . , r„. The later symme- 
try leads also to certain type of symmetry for the permutations of momenta pi, . . . ,Pn, as it 
follows from Bose-properties for whole amplitude M. Notice that the charge distribution differs 
from Poisson one if amplitude M is symmetric relative to permutation of momenta. As it will 
be seen below, we get in this case the significant increase of the probability of events with 
neither neutral nor charged pions. On the other hand, the permutation properties of amplitude 
M can give very useful information on properties of the objects generating the Centauro events. 

For a system of n pions with small values of total isospin the charge distributions for 
certain types of symmetries (Young tableaux) were derived more that twenty years ago by 
Pais The charge distribution in the case, when all possible Young tableaux give the same 
contribution to the wave function, was obtained in [0. In this case the distribution has sharp 
form and is close to Poisson one for large values of n. 

^In ref. |Q] concrete realizations for all Young tableaux with one and two rows were constructed, some special 
cases with three lines were considered also. In connection with this question see also paper pi. 



1_ 

no = -n 



2_ 

rich = -n- 
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Recently in connection with Centauro events in ref. [0 the role of small isospins in a 
system of n pions was emphasized. This is an important fact because in the nuclear collisions 
namely small vlues of isospin are realized. To explain Centauro events in ref. [|^ the most 
symmetric form of wave function for a system of n pions was also used. These two assumptions 
automatically lead to the smooth charge distribution, in accordance with ref. 

In the present paper we are going to collect altogether those results of refs. [§, 0] , which 
in our opinion may take immediate attitude toward Centauro problem. Transparent deduction 
of the results will be given. Taken altogether, these results give rigid restrictions on 

theoretical models of Centauros. We consider this paper as a short review. 

Physical reasons will be explained, which enhance the production probability of the pions 
of the same sort. This enhancement takes place when the totally symmetric (in momentum 
space) amplitude dominates, and it compensates Poisson's diminition of process. This phe- 
nomenon has much in common with induced photon emission, that leads to formation of laser 
bunch. 



2 Connection between the symmetry of amplitude and 
the charge distribution of pions 

The production amplitude of n pions has the form: 

M{pi,Ti;p2,T2;... ;p„,r„) = ^Mo(pi,... ,p„)v9^(ri, . . . ,r„), (1) 

a 

where the sum is taken over all the possible types of symmetry a, characterized by Young 
tableaux. As we have only three sorts of pions (vr^, vr", vr"), index {i = 1,2, ...,n) takes also 
only three nonzero values. It follows that the rows number in Young tableaux doesn't exceed 
three as well. Totally symmetric function ipa^Ti, . . . , Tn) corresponds to the only Young tableau 
with one row. 

The most important fact consists in the following: more symmetric function (pa{Ti, ■ ■ ■ , Tn) 
gives more smooth distribution in charge space. On the contrary, less symmetric Young tableaux 
with three rows give sharp charge distribution. Notice that the total number of young tableaux 
is very large: 

3 /s" 3*^ 



and most part of comes from Young tableaux with three rows. So to observe the Centauro 
events with appreciable probability, it is necessary to increase the role of symmetric terms (with 
one and two rows) in the sum (|1]). In other words, it means that by some dynamical reasons 
the pions are produced in symmetric states in momentum space. 

First of all, let us demonstrate the pion distribution in the case when large number of 
amplitudes with the different symmetry types gives the same contribution into the sum (|l|). 
Consider the case of even number n = 2k with total zero isospin / = 0. We shall obtain the 
pions distribution for the amplitude in a special form: 

M = S {M{pi,p2, . . . ,Pn){ava2) ■ ■ ■ {a2k-va2k)} , (2) 
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where Si (i = 1, . . . ,n) is isotopic vector of pion, S is the symmetrization operator. Equation 
(§) for M is the special form of (|l|), as eq. (H) is symmetric relative to permutations of pions 
in pairs, e.g., ai 02, etc. Nevertheless this form of M gives rather sharp distribution in 
number of neutral or charged pions. Suppose that M{pi,p2, . . . ,Pn) depends sharply on the 
arguments pi,P2, ... ,p„ so as it differs noticeably from zero only when pi ^ Pi,P2 ^ P2, ■ ■ ■ 
and Pi ^ P2 p3 ■ ■ -Pn- This condition provides that the amplitude M, eq. (|), contains 
a large number of terms with different permutation symmetries. Making symmetrization, we 
get M{p2,pi,P3, ... ^ for pi ^ pi, p2 ^ p2- On the contrary, M{p2,pi,P3, ... ,p„) differs 
from zero when pi ~ p2, P2 ~ Pi, Ps ~ Ps, • • • • We get the same result making permutation of 
any pair of arguments. So, integrating over whole phase volume to calculate the cross section, 
we obtain that all the permuted terms in eq. (H) give the same contribution. It follows that 
the value of cross section is proportional to a number of nonzero isospin amplitudes in eq. (0). 

Now let us find out how this number depends on the number of charged (or neutral) 
pions. Assume that among n = 2k pions there are ki of 7r+, ki of vr" and 2^2 of vr^-mesons. 
Considering vectors in (0) in the coordinates a±, Oq, we get for the isospin part of amplitude: 



X 



A-n ... A " 



X-^i ... A " 



(3) 



As it is seen, ^ = 1 for ri = . . . = = 1, r{ = . . . = r^^ = -1, cxi = . . . = cr^a = . . . = = 0- 
Let's look for such permutations of mesons, which give us nonzero values of x- First we can do 
all the possible permutations of upper indices in the first brackets in (H), the set of lower indices 
being fixed, and obtain the factor ki\. Secondly it is also possible to permute indices in second 
bracket in eq. (|^). This operation changes positions of 7r°- mesons (no identical operations have 
been done!). The number of nonzero amplitudes is {2k2 — 1)!!. Permutations of indices between 
both brackets give zero result.^ From this it follows that the probability w{ki,k2) to emit ki 
pairs of vr+vr^ and k2 pairs of 7r° is 

The denominator in (^ appeared from the usual normalization condition for the cross section 
with an identical particles in final state, see, for instance, ref. Calculating the normalization 
constant for w{ki, ^2), we finally get: 

k\ f2V' fl^^^ 



where k = ki + k2- This is binomial distribution. For large values of k and small ki (or ^2) 
we may approximate eq. (^ by Poisson distribution. It follows from eq. (^ that the emission 
probability of nothing but charged particles only (i.e., the probability to observe the Centauro 
event) is suppressed strongly, e.g., ^(Tr+Tr") ^ 10~^ for n = 100. 

In the deduction given above we used for amplitude M partially symmetric form, eq. (|^). 
Surely there exist other combinations of isovectors Oj of less symmetric form, e.g. ay [a2 x a^]. 

^Total number of terms in (||) and in the cross section also is 2*''fc!A:i!(2fc2 — 1)!!. Later on we shall study the 
distribution over fci and ^2 only. 
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Being included in M, they would give more sharp charge distribution then it follows from 
In the most consistent form tetrms of all different symmetries were taken into consideration 
in ref. 0. Contribution to the cross section from all that terms was taken to be equal. In 
comparison with (|^) the asymptotics of formula (16) from ref. gives the probability of event 
without charged pions 



Analogously obtained probability of nothing but charged particles has the form: 

This is less than it follows from eq. (^) at n = Ik. 

Let us consider now the contrary case. We find the charge distribution in a system of 
n = 2k pions at J = 0, described by the state totally symmetric in the momenta (and isospin) 
variables. Corresponding isospin wave function is obtained from eq. (|^) provided the amplitude 
M{pi, . . . ,pn) is not changed under the momenta permutations. We emphasize that in this 
case the amplitude (|^) is not yet the amplitude of a particular form. In this case it is defined 
unambiguously and doesn't depend on the angular momentum addition scheme in the initial 
amplitude which produces after symmetrization the amplitude eq. (0). Therefore all nonzero 
terms in eq. (§) are the same, that results in 100% constructive interference. Then instead of 
eq. (^) we obtain 

Calculating the normalization factor by means of the relation: 

(2Z-1)!! (2A;+1)!! 



(a)!! (2k)» ' 

where, by definition, (—1)!! = 1, we find: 

, X ^''kl (2A;2-1)!! 

In contrast to eq. (^, the distribution (|^) is rather smooth. It is shown in figure [I|. At 
/c2 ^ 1 we have w{k2) ~ l/v^- The average numbers of the neutral and charged pions both 
in the case of distribution (|^) and in the case of eq. (^ are found the same and equal to 

1 2 

Uq = -n, rich = -n. 

/^From eq. (0) at n ^ 1 we obtain the production probabilities of nothing but neutral pions 
w{tt^) or nothing but charged pions w(7r+7r^) (Centauros): 
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Figure 1: The function y(m) = (2m — l)!!/(2"*m!). It gives multiplicity distribution over a 
number of paires of neutral pions. 



These probabilities are enhanced by many orders in comparison to eq. (|^). This enhance- 
ment is just well known in quantum optics the induced emission phenomenon. The numerator 
of eq. d^) contains extra (in comparison with eq. (^) factorials ki\ {2k2 — 1)!!. These extra 
factorials appear due to interference, which is extremely essential in the case of symmetric state 
relative to the final pion momenta permutations. The same behaviour of amplitude is found in 
the case of pion emission with the same momenta. The reason of this enhancement results from 
the fact that the emission of one sort bosons in the state with the same momenta has essential 
advantage over the emission of the same number of bosons of two or several sorts. Note that 
according to eq. (^) the probability w{ti^) by (7r?7./2)^/^ times less than w{n'^n~), in contrast 
to the Poisson distribution. In the latter case we have: 

However, inspite of maximal enhancement of neutral pion production probability ^(vr*^), it 
remains less than w(7r+7r~). This fact is connected with particular diminution of "initial" 
Poisson probability w{tt^). 

This can be also explained without starting from the Poisson distribution. Although the 
enhancement of induced vr^-emission is considerably stronger than that in the case of the same 
number of tt+tt" (see below eq. (p!OD), however the combinatorial factor (the coefficient 2) 
enhances the vr+vr" -pro duct ion (see below eq. (|ll|)). This factor appears in expansion of the 
scalar operator c^^ = 2c^ + c^_ + c^o in terms of the pion creation operators. 
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In order to emphasize more clearly the connection between possible Centauro production 
and the induced photon emission we reproduce eq. (|^) by means of the second quantization 
formalism. The totally symmetric isospin n-particle wave function at / = has the form: 

(c^^)'lO), 



where c^^ = 2c|^+c|^_ + c^^o, c] is the pion creation operator of i-th sort, n = 2k. Calculating 
the amplitude (kiir'^ , kiiT~ ,2k2TT^\(f)) and taking into account that 



{k,n^\icU)'^\0) = Vh^. , (2A:27r°|(cioriO) = v/(2^, (10) 

we find: 



(M+,A;ivr-,2A;27r'^|0) = Cf ^^=^A^, (11) 

V(2A; + 1)! 

where C^^ is the binomial coefficient. The amplitude (|TTp squared reproduces eq. (|]). En- 
hancement of the pion production amplitude in the case of the same sort of pions is due to the 
same factors ([To|) which enhance the photon induced emission amplitude. 

Note that the formulae (^ for w{it^) and w{7i^tt~) coincide with the result obtained in 
ref. 03 with accuracy of replacement of the pion number n in eq. (D by the average number n 
in ref. [0. The probability of a given number of pions in ref. [0 was described by the Poisson 
distribution. Averaging eqs. by Poisson distribution just leads to the replacement of n by 
n (with accuracy of the power corrections). 



3 Role of the small total isospin of the n pion system 

It was shown in the previous section that refusal from domination of most symmetric isospin 
wave function leads to negligibly small probability of Centauro production. In the present 
section we show that refusal from small isospin values also sharply diminishes the probabili- 
ties ^(vr'') and w{tt'^tt~) even in the case of totally symmetric wave function. This aspect of 
phenomenon has been emphasized in ref. 0. 

Les us calculate the probabilities w{7r^) and w('7r"'"7r~) considering the totally symmetric 
n-pion state with n = 2k and isospin I = n, I3 = 0.0 The isospin wave function has the form: 

|0,) = A(/)(ct^)'=F,o(-^)|O), (12) 



where ^4(7) is the normalization factor: 

47r 



A{I) 



(2/+1)!!' 



■^In ref. [Q the probabilities w(7r*') and w{tt^tt^) were calculated by projecting pion field on the state of 
small isospin, / = 0, 1. 

''The charge structure of the state with J3 close to the maximal value is almost definite. Therefore the 
problem appears for the state with small I3 only. 
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The probabilities concerned are determined by the following amplitudes: 



M(7r°) = (2fc7r'^|0,) = A(/)y^^P,(l)(2fcvr°|(c!;o)^1O), 



where Pi{x) is the Legendre polynom. 
^From here at / ^ 1 we obtain: 



w{7r') = (13a) 
^(vr+vr-) = ^. (13b) 



At / = n ~ 100 the probabilities (|13a|) , ( |13b|) have the order of 10 



1-30 



4 Conclusion 

Thus we have convinced that the events containing nothing but 7r"'"7r~-mesons or vr'^-mesons 
at high multiplicity can have appreciable probability only at simultaneous realization of the 
following two conditions: 

1. Symmetric state of n-meson system relative to the permutation of momenta (and isospin 
indices) . 

2. Small total isospin. 

The first condition seems us to be considerably more essential than the second one, since 
due to isospin conservation the second condition is satisfied automatically, while the first con- 
dition implies the strong restrictions to the dynamics of high multiplicity pion production. We 
emphasize that these results follow from rather general consideration. Therefore any model, 
which pretends to explanation of Centauros, must first of all satisfy to two these conditions.^ 

It should mention that under symmetric function we mean the function corresponding to 
the Young tableau with one or two rows. Thus, proceeding from ref. 0] one can obtain that 
/ = for the Young tableau with two equal rows w(7r+7r^) = 1/n, but w(7r°) = 0. However, the 
Young tableau with three approximately equal rows leads to distribution which differs from zero 
only at uq and rich close to |n and |n correspondingly. Since the overwhelming majority of the 
Young tableaux is the Young tableaux with three rows, namely they determine the distribution 
when all the Young tableaux give equal contributions. 

On these grounds it can be conceived that Centauros are produced from decay of a particle 
as if consisting from large number of pions in totally symmetric state. On the ground that the 
experimental hadron transverse momenta in Centauro equal to pt > 1-5 GeV/c, one can expect 
that the size of this hypothetical object is of the order of (1.5 GeV/c)~^. The field with energy 
E ^ 230 GeV is concentrated in this volume. Considering the pions as quanta of this field 

^We mean the models in which the hadrons produced are considered as pions. 
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one can suppose that this object could be similar in its properties to the classical solution of 
a field equation. One can expect that namely at high energy the heavy particles having the 
properties close to the classical solutions can be produced. This hypothesis have been suggested 
in ref. 0. Note that possible soliton solution to explain Centauros is, for example, the slowly 
damping classical solution for Higgs scalar field found in ref. [0. 

Let us give a possible estimation of the Centauro production cross section. Assuming 
that a Centauro is produced from a soliton decay, we see that the only dimensional parameter 
is its size ~ (1.5 GeV/c)~^. It follows that relative probability of the soliton production is 
the ratio: 

w ~ rl/rl ~ 10-^ (14) 

where Vst is the strong interaction radius. This estimation is rather rough. The probability of 
Centauro production is obtained by multiplying eq. ([I^) by (^. Therefore the probability to 
find a Centauro is estimated as 

u;(Centauro) ~ 10"^ (15) 

Although this estimation does not depend explicitly on energy, it becomes valid only 
begining from rather high energy (which we don't estimate), when quasiclassical objects can be 
produced. CERN SPS-coUider experiments |TI|, [T^ at ^/s = 540 GeV give so far no evidences 



for Centauro production. From the present point of view this means that inspite the fact 
that the energy ^/s = 540 GeV exceeds the threshold for production of the mass M = 230 
GeV/c^, this energy is not enough for proceeding the physical phenomena leading to Centauro 
production (e.g., for developing the soliton object). 

The other approaches to Centauros were discussed in refs. 0, O . 



5 Appendix 

From the present consideration some conclusions follow for the high multiplicity pion production 
in nuclear reactions. It follows that charge distribution of pions produced in any nuclear reaction 
would be smooth provided the pion momenta are close enough to each other. The latter comes 
true in kinematics near the phase volume bound. For example, in the reaction 

pp ^ pp + rnr (16) 

at highest possible pion number n allowed by energy conservation the reaction kinematics 
automatically forbids the wide variation of the pion momenta. This ensures automatically 
the total symmetry of pion production amplitude relative to the momenta permutation in the 
narrow allowed region. As it was shown above, we obtain the smooth charge distribution in 
this case. 

It would be useful to carry out the experimental research of the reaction (0) in the 
kinematics concerned. Such investigation is possible at the meson factory accelerators, where 
due to high intensity the small cross sections can be measured. 

Another interesting reaction is the extremely many pion annihilation of low energy an- 
tiprotons (see ref. 0]): 

p + p nn. (17) 
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Usually rather sharp binomial distribution (as eq. (§)) is expected in such processes, 



which follows from the statistical models |]T5l. In this reaction at highest possible multiplicity 
n we expect a smooth distribution as well. 

The authors are grateful to A.B. Kaidalov, V.E. Markushin, A.I. Nikishov, V.A. Novikov, 
L.B. Okun, I.S. Shapiro and M.B. Voloshin for useful discussions and valuable remarks. 
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